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Abstract. We demonstrate the transport of interface states in the one- 
dimensional ferromagnetic Heisenberg model by a time dependent magnetic field. 
Our analysis is based on the standard Adiabatic Theorem. This is supplemented 
by a numerical analysis via the recently developed time dependent DMRG 
method, where we calculate the adiabatic constant as a function of the strength 
of the magnetic field and the anisotropy of the interaction. 



1. Introduction 

The problem of calculating the dynamics of quantum spin chains has acquired 
new importance due to potential new applications in spintronics [1] [2], quantum 
information, computation and control theory and the now realistic possibility of doing 
experiments on systems that are accurately described by a one-dimensional array of 
spins [HI 11] . The development of time-dependent versions of the successful Density 
Matrix Renormalization Group (DMRG) method of White [HI H] by several authors 
[3 El m [12 therefore comes at a fortuitous moment. All numerical schemes for 
a quantum many-body problem such as a spin chain necessarily involve a drastic 
reduction of the high-dimensional Hilbert space to a suitable subspace of relatively 
modest dimension. The DMRG method does this by approximating the desired states 
(typically the ground state and low-lying excitations above it) by finitely correlated 
states also known as generalized valence bond states or matrix product states. 
The time-dependent DMRG method extends this idea by allowing the subspace used 
in the approximation scheme to depend on time. 

Here we are interested in the dynamics of the Heisenberg XXZ chain with kink 
boundary conditions [12j and a time-dependent external field. The ground states of 
this model are the so-called kink states which describe an interface between domains of 
opposite magnetization. Their degeneracy grows with the length of the chain [T3lll4j . 
and is labelled by the position of the kink, or equivalently, the total magnetization 
which is a conserved quantity. A transverse external field localized at one site can be 
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used to select a unique ground state |15| . i.e., such a perturbation pins the kink at 
a specific site. These interface structures also appear in certain asymmetric simple 
exclusion models of stochastic particle dynamics [Tni[T71[IH] and it is a natural question 
to study their non-equilibrium properties. In this paper we study the quantum 
dynamics of the kink under the influence of a time-dependent magnetic field. Typically 
we use a field supported at one or two sites and moving at a constant speed. For small 
velocities, we know by the adiabatic theorem that the time-evolution of a ground 
state will follow the ground state of the time-dependent Hamiltonian. The standard 
heuristic scale to identify the adiabatic regime is the smallness of the ratio of the 
time-derivative of the Hamiltonian over the spectral gap squared. We will verify the 
usefulness of this criterion by calculating the time evolution of kink states in the XXZ 
chain under the influence of a time-dependent magnetic field. Needless to say, the 
dynamics of kink states, which can be regarded as discrete analogues of solitons, is of 
great interest in its own right. Moreover, they provide a simple dynamical quantum 
model of sharp magnetic domain walls. 

In order to numerically study the quality of the adiabatic approximation and its 
domain of validity, it is necessary to express the evolved state in the time-dependent 
eigenbasis of the Hamiltonian. Therefore we construct the DMRG basis at time t 
starting from the ground state(s) and low-lying excited states of the Hamiltonian of 
the system at time t. This is in contrast with the original method, where the time- 
dependent DMRG basis is obtained from the evolved state at time t. Our method 
can be of more general interest, since the calculation of the DMRG basis is a separate 
calculation about which we can be confident to have good control, regardless of the 
potential difficulties in reliably calculating the time evolution. 

The paper is organized as follows. In Section [21 wc define the model and study 
its basic properties. The adiabatic approximation is studied in Section [3l In Section 
|4j we give a brief discussion of the issues that arise with rapidly changing magnetic 
fields. Details of the DMRG algorithm are given in [Appendix A| 

2. The model 

The spin-i ferromagnetic XXZ Heisenberg model on the chain [1,L] with interface 
boundary conditions is defined by the Hamiltonian Hq ~ X)f=i hx,x+i, with nearest- 
neighbor interaction 

^ -^-\SlSl+, + SlSl^^) - SlSl^^ + yi-/^-^{Sl ~ Sl_,^) + il , (1) 

where A > 1 is the anisotropy parameter and the matrices (S*^, 5^, S^) are the usual 
spin-i matrices; hj.^j^x is a projection and has ground state energy zero. We have 
set h = 1 which should be kept in mind when we define adiabatic regimes in terms of 
the magnitude of the magnetic field and the velocity. 

The Hamiltonian iJo acts on the Hilbert space TLl ~ 'CP' ■ It has a large 
(quantum) symmetry with a similar multiplet structure as the SU (2)-symmetric, 
isotropic model where A = 1. In the limit A — > oo we obtain the Ising Model, 
Rising _ i/o(A = oo). The ground states and excited states (for finite L and 
in the infinite volume limit L ^ oo) have been analyzed completely in recent 
years [HI [2D1 [HI [521 US]- A convenient basis of ground states are the so-called 
grand canonical ground states, (pc- They are product states and depend on a complex 
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parameter c G C \ {0}: 



(1 + icw^^r^' ( \. ) 



(2) 

The state (jjc describes an interface state that is exponentially localized at xq — 
— ln|c|/ln(7 with q ~ A — VA^ — 1 G (0,1). The width of the interface depends 
on A and becomes sharp (i.e., the transition from up to down spin occurs across one 
bond) in the Ising limit A ^ oo and flat as A | 1. There are no interface states in 
the isotropic Heisenberg model. 

We perturb the Hamiltonian by a magnetic field B : [1,L] — > M'^ satisfying the 
locality condition: 

(Al) The support of B is finite, non-empty, and independent of L. 
Then we define the Hamiltonian 

Hv =Ha + V with V = Y^ B{x) ■ . (3) 

X 

We refrain from considering the infinite chain limit in detail. Condition (Al) ensures 
that Hv is still a well-defined semi-bounded operator in the GNS representation of 
the unperturbed Hamiltonian. 

The case of a magnetic field located at a single site y has been analyzed in [TS] . 
Let B1+1B2 7^ Q,Cy = -{\\B\\+B3){Bi+iB2)~^q~y, and let be a state of the form 
([2|). This state (pc^ has energy — i||Bj|, and since the ground state energy can be at 
most shifted by this amount, we know that it is a ground state. It is also the unique 
ground state which moreover has a uniform gap as i — > 00. If B = (i?,0,0) with 
B 0, then (pcy is a kink localized at site y. If B(a;) = [0,0, B{x)) with B{x) < 0, 
then the ground state is the all spin-up state which, as L — > c» is however not in the 
GNS kink Hilbert space; see the discussion in [15]. 

In the Ising limit A 00, ground states may be largely degenerate as for V = 0. 
For example, consider the magnetic field B(l) = (1,0,0),B(2) = (0,0,1), and zero 
otherwise. Then Hv has L — 1 "kink" ground states of the form '0[i,2] ® ipis.L]- 
Here, '0[i.2] is the ground state of the two-site Hamiltonian which has the form 
V'[i,2] = (0,c,0,-Vl - c2) for some c > 0. Note that (V'[i,2] |5'2'0[i,2]) = -1/2- Thus 
we can choose V'[3,l] to be an (Ising) kink state on the remaining sites of which there 
are L — 1. 

Thus we consider fields B such that for all x with B(.t) 7^ 0, B(a;) has a non- 
vanishing component in the xy-planc. This means that, in addition to (Al), we 
assume the following: 

(A2) On the support of F, B(.t) 7^ (0, 0, ±||B(.t)||). 

Theorem 2.1. Let B satisfy the conditions (Al) and (A2). Then, there exists a 
finite Ao = Ao(B) so that for all A > Aq, 

(i) Hv defined in (0) has a unique ground state; 

(ii) Hv has a positive gap above the ground state uniformly in L. I.e., there exists a 
7(A) > 0, independent of L, such that for all states ip orthogonal to the ground 
state we have (V'KiJy — i?o)|^) > 1 II^IP- 
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The conclusions of the theorem can be expected to be vahd for other situations 
than the ones covered by (Al) and (A2). In the cases where we have performed 
numerical calculations, a positive gap appeared for all A > 1. Here we only 
show the result for sufficiently large A by applying the theory of relatively bounded 
perturbations, see e.g. [24]. Let us recall that an operator B is relatively bounded with 
respect to the operator A if there exists a (finite) constant M so that for all in the 
domain of A, \\B^\\ < M(\\Aip\\ + Then, all eigenvalues and eigenprojections 

in the discrete spectrum (i.e., isolated and finite multiplicity) of A{a) = A + aB are 
analytic for a £ {—ao,ao) for some ao > (cf. [Ml 3.5.14]). 

Proof. Let Hy"^^ = Hv{A = oo) be the Ising kink Hamiltonian with magnetic field B. 
By using a rotation around the 2;-axis we can assume that 'B{x) ~ {B^{x),0, B^{x)) 
with B^{x) < for all x in the support [a,b] of B. Then the Hamiltonian Hy"^^^ 
defined on the support of B has a unique ground state '0[a,b] by the usual Perron- 
Frobenius argument. Le., there is an rt G N such that (^^y'["\])" has strictly negative 

off-diagonal entries. We claim now that = '^x=i (i) ^ V'[a,b] ® ^x=b+i (o) 

the ground state of iJy™^ on the full chain [1,L]. By adding the other (non-negative) 
Ising Hamiltonians to the left and right of the support of B the ground state energy 
could in principle increase. But fpii.L] is also an eigenvector of i/^'"^ with the same 
eigenvalue as "0(0, &] since the additional terms vanish on Thus ^[i.l] is the 

ground state of iJy'"^. 

Next, we consider the Heisenberg Hamiltonian iJy(A) = if^™^ + A~^K -f (1 — 
— A~^)P. P only involves the boundary spins and is therefore uniformly bounded. 
K contains the XY terms in the Hamiltonian and the norm of this term is of order 
L. In particular, it is not uniformly bounded but one can show that K is relatively 
bounded with respect to iJy'"^. This was done in |5S] for the Hamiltonian without 
i?-field. The relative bound is uniform in L. By our assumptions, the i3-field defines 
a bounded operator and hence K is also relatively bounded with respect to Hy"^^. 
It follows that for A~^ in a non-empty interval [0, A^^], the spectral projection for 
an interval of energies containing the unperturbed ground state energy is an analytic 
function of A~^. This implies the statement of the theorem. □ 

Our prime example are magnetic fields located on two adjacent sites, xq ^ L/2 
and xq -I- 1. More precisely, let u, i? > and define 

B{x,t) ^ -B ■ f{x-vt) , 

)l-|.T-a:o| for cco - 1 < X < cco + 1 

"In ■ 
I otherwise 

The perturbation 

L 

V{t) = ^ B{x, t)Sl = -{l-vt)B Sl„ -vtB Sl+, (4) 

x=l 

for < t < T ^ 1/v satisfies the conditions (Al) and (A2). Hence, we know from 
Theorem 12. II that the Hamiltonian, 

H{t) = Hn + V{t) (5) 

has a unique ground state, ip{t), with a positive gap. The ground state energy 
Eo{t, B) is analytic on (0, v~^), also as L ^ oo. It is well-known [24l 3.5.23] that the 
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ground state energy EQ{t,B) is concave in t. By symmetry we also have i?o(i, -B) = 
Eo{v-^ - i, S) for < < < l/{2v). This implies £^0(1/2^, B) = maxj Eo{t, B). We can 
get more information on EQ{t,B) from the low energy spectrum of the Hamiltonian 
with the magnetic field at a single site. For simplicity we put v = 1. 

Proposition 2.2. For t > let g{t) be the gap above the ground state (p^y of 
Hy(t) ~ Hq — tBSy. Then g{t) is increasing in t (strictly speaking, g{t) also depends 
on L which we will tacitly ignore). Further, on [0, 1/2] we have for the Hamiltonian 
H{t) defined through ^ and ([^, 

i?W>^5(2t)(l-|(0cJ0c„+,)|). 

Proof. We may assume that i? > 0; otherwise apply a rotation. That g{t) is increasing 
follows from the variational principle. To this end, \eiO <t <t' . Let ip be orthonormal 
to 4>cy such that {ip\ Hy{t') + \t' B = g{t'). Then, 

(^1 Hy{t) + Ub IV.) = g{t') + {t- t')B (V-l \-Sl IV') < g{t') . 

Thus, g{t) < g{t'). The estimate on the ground state energy of H{t) is based on 
Kitaev's Lemma [26\ Lemma 14.4]. In general, let Ai and A2 be two operators on 
some Hilbert space Ti. By Gi_2 and Hi. 2 we denote the space of ground states and the 
orthogonal projections onto Gi,2, respectively. We assume that Gi n G2 = {0}, and 
that for some S > 0, Ai^2 > S{1 - 01,2). Then, 

Ai+A2>5(l- sup |(Vi|V'2)|). 

For our proof it suffices to apply this inequality with 

Ai = ii/o - (1 - t)B{Sl - il) = iffj,(2(l - t)) + 1(1 - t)Bl > , 
A2 = \Ha - tB{Sl+, - il) = ^Hy+,{2t) + \tBt > . 

□ 

We supplement this analysis of the ground state properties of the family of 
Hamiltonians ([5]) with numerical DMRG calculations for a chain of length i = 20 
and an anisotropy of A = 2. The chain length is chosen such that for given A, L is 
much larger than the width of the interface, but otherwise there are no qualitative 
differences for different combinations of A and L. We distinguish three ranges of 
B- values with qualitatively different ground state properties. 

A typical intermediate value of B is 0.5. The 5''^-magnetization profile of the 
ground states Lp{t) of H{t) nicely interpolates between the interface product states 
(p{0) = (f)^^^^ and Lp{l) = (/>c,„+i (see eq. ^) with .tq = L/2 (Figure [1]). The first 
excited states are still well localized around the interface positions. Both the groimd 
state and first excited state energies are concave as a function of t (Figure [5]). The 
ground state energy is minimal and equal to the known value — _B/2 at the start {t = 0) 
and end {t — r) points. The gap is always strictly positive, but minimal at half period 
(t = T 12). The transition from 1^9(0) to 1^9(1) is not homogeneous, in the sense that the 
right or top half of the interface moves first and the left or lower half moves later. This 
effect becomes more pronounced as B increases (compare Figure [T] left and right) and 
can be understood as follows. In the first phase after switching on a magnetic field at 
xi^ + 1, »S'^o+i has the effect of rotating the up vector at xq + 1 into the (1, 3)-plane. 
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Figure 1. S'^-magnetization profiles of the ground states for B = 0.5 (left) and 
B = 5 (right), for < t < r/2 (■) and t/2 < t < r (• ) (r = 200, 200 time steps, 
profiles plotted for every 20 time steps, only central sites around the interface are 
shown). 




Figure 2. Energy of the ground states and first excited states (top) and energy 
gap (bottom) for B = 0.05, B = 0.5 and B = 5 (left to right). 

But here nothing happens at site xq because the initial state is an eigenstate of S'j.^^ . 
The same thing happens at the end of the cycle when the vector at xo + 1 is almost 
an eigenvector of •S'^o+i- Site xo is affected only after some time as the perturbation 
at xo + 1 is communicated to xq in second order. 

For very small magnetic fields, the ground state energy is still a concave function 
of t, but the first excited state energy is now convex (Figurc[2]). The energy gap is still 
minimal at half period. Furthermore the energy gap for a 10 times smaller magnetic 
field (0.05 vs. 0.5) is approximately 10 times smaller, confirming the theoretical result 
[T5] (for a single site impurity) that the gap scales linearly with B for small B. For 
magnetic fields much larger than the intermediate value B ^ 0.5, the phenomenon 
that the interface moves in separate steps becomes much more pronounced (Figure[T|). 
Like for intermediate i3- values, the energies of the ground and first excited states are 
concave functions of t, but the minimal gap now occurs at the start and final times 
(Figure El). 

3. Adiabatic transport 

Let us recall the standard adiabatic theorem which applies also to the case of an 
infinite chain. Let r > and H{t), < t < t, he a. family of self-adjoint operators 
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with common dense domain, and let tp{t) be the solution to 

= rH{t)m (6) 

with initial condition iJj{0). By P{t) we denote the spectral projection onto the ground 
states of H{t). We assume that P{t) is piece-wise, twice continuously diffcrentiable, 
finite dimensional, and uniformly (in t) separated from the rest of the spectrum of 
H{t) by a gap 7(t). 

Theorem 3.1 (Adiabatic Theorem, Kato). Under the above conditions on H(t) and 
ip{t), there is an eigenvector ip{t) of Hit) with ip{Q) = ipifi) and a constant C such 
that, 

sup U{t)-^{m<CT-\ (7) 

0<t<r-i 

Proof See [57] or [H 3.3.11]. □ 

We call the smallest constant C in ([7]) the adiabatic constant. Heuristically [HI 
17.112], the adiabatic constant is of the order 

^^P rvz — ' S^t T . (8) 

sG[0,l] l{.sY 

In our case, H{t) is given by (P, and \\f^H{s)\\ = \\BSl^^ - BS'^„+i|| = B. Hence, 
C (X i?/7^i„ with 7niin = minsg[o 1] 7(s) is the minimal gap. As B tends to infinity 
the gap 7 saturates and therefore C grows linearly with B for large B. On the other 
hand, if B is small then the gap shrinks like B (not B^) as has been shown for a single 
site perturbation [15] . and C diverges as l/B. As a consequence, there is an optimal 
range for which C is smallest. In Figure[3]this appears near B ~ 1. If we fix B and A 
and some e we can find (empirically) an upper bound Wmax for the velocity v for which 
then the adiabatic evolution is e-close (in the £^-sense) to the true time evolution. The 
slowdown (decreased Wmax) of the domain wall motion for large B-fields was reported 
cf. [2]) but here we also observe a slowdown for small S- fields. On a rigorous level, 
estimates in the vein of ([8|) on the adiabatic constant were recently derived by Jansen, 
Ruskai, and Seller 

We have computed the adiabatic constant numerically using adaptive time- 
dependent DMRG [71 in] . Unlike the original method, we express the evolved state 
in the time-dependent eigenbasis of the Hamiltonian such that we can compute eq. 
(jT)) with ij^it) a-nd (p{t) expressed in the same basis (see [Appendix A| for algorithm 
details). Figure [3] shows the adiabatic constant as well as the heuristic upper bound 

B 

C < a— , 

with the energy gaps computed by ground state DMRG (sec Section [2]). To determine 
a, we took one value of B {B = 1.0 for A = 2 and B = 2.0 for A — 10) and made the 
inequality in ([5]) an equality at the lowest v {v = 0.001) at this particular value of B. 

Since our time-dependent algorithm does not target the state ^(0 directly, but 
rather the lowest energy states of the Hamiltonians H{t) (see [Appendix A[ ), it is 
important to keep track of how well ip{t) is represented in these time-evolving DMRG 
bases. One way to measure this is by computing the deviation from 1 of the norm of 
ipit). As expected there is a correlation between this norm loss and the value of the 
adiabatic constant C. In Figure [3] where C is minimal, ||V'(OII '^i^l be above 0.999 at 
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Figure 3. Adiabatic constant as a function of magnetic field strength B for 

diff'ercnt velocities v = 0.01,0.005,0.002,0.001 ( , bottom to top) compared 

to the upper bound c*S/7^in ( ) for anisotropy A = 2 (left, a = 0.5677) and 

A = 10 (right, a = 0.6324) with L = 20. 



all times; at the very smallest and very largest i3-fields, and for the largest velocities 
still considered adiabatic {v « 0.01) the minimmxi norm drops to about 0.97. 

Adiabaticity is also measured by comparing the value of the total energy in the 
time-evolved state with the ground state energy of H{t), i.e., 

should be close to at all times. For a typical adiabatic speed v = 0.005 and 
intermediate -B-field, this difference never exceeds 5 • 10^"'. Energy differences in 
the small and large B regions are of similar magnitude. 

Finally, we comment on the role of the interpolation function ^ in our adiabatic 
Hamiltonian. Let us consider a large magnetic field. Then as time proceeds part of 
the profile is lagging behind (cfr. Figure [T]). Also, the change of the profile is rather 
quick in the beginning (and at the end) and slow in the middle of a cycle. To improve 
on the transport of the domain wall for a fixed time interval we can take advantage of 
this phenomenon and slow down the time evolution in the beginning and accelerate in 
the middle by choosing a different interpolation function. In other words, we may use 
a general interpolating Hamiltonian Hf{t) ~ H^ — {1 — f{vt))BS}^^ — f{vt)BS}^^_^_i with 
the constraint that /(I) = 1-/(0) = 1 to keep the same mean velocity. As an example 
we have used the function fi{t) = cos(7rt/2) and a piece- wise linear function /2 with 
slope 1/3 for t e [0,1/6] U [5/6,1] and slope 4/3 on the interval [1/6,5/6]. In the 
first case we took L = 20, A = 2, B = 5,v = 0.05 and reduced the adiabatic constant 
from 16.6 (constant velocity, or equivalently linear /) to 4.5 calculated with /i. In the 
second example with /2 and parameter values i = 20, A = 2,_B = 10,1; = 0.05 the 
adiabatic constant dropped from 25.1 to 12.2. This observation is in agreement with 
|29] . where the adiabatic constant for the general interpolating Hamiltonian Hf[t) was 
studied. 



4. Fast change of the magnetic field 

Now we study the situation when the velocity of the moving magnetic field is large. 
We start initially in the ground state cj) = (j)c^^ of the Hamiltonian H{Q) — Hq — BS].^. 
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Let T = 1/v and let p(t) be the probability that at time r the system is still in the 
state (/). Then according to formula [SS] (17.60)], 

pir) = l-T'var^(iJ)+0(T3) 
with va.T^{H) = {(t)\H^\(j}) - {(l)\H\(l))'^ , and H = ^ H{t)dt. In our example, 
H = Hq — -f ("^xo + ^xo+i)- Since (j) is also an eigenstate of Hq — -f 5^^, with energy 
— we obtain 

Hence, 

ya.T4H) = ^{l-A{q^\Sl^^,\q^f). 

The quantity {<j>\Slg_^_l\(j)) can be computed fairly explicitly as a function of A but 
we are content with the trivial estimate that 1 — 4((/)|S'^jj^j^ |(/))^ < 1. As a result, the 
probability to stay in the initial state until r, 

p{r)>l-^ + 0{T'). 

If we want this to be larger than 1 — e, then \B\ needs to be smaller than Av^/e. 

Numerically we can also investigate the intermediate region between adiabatic 
and sudden change of the magnetic field. Here, B and v are of the same order of 
magnitude. In Figure IH we follow the S'^-magnetization profile of the time evolved 
state over two periods, with the natural extension of V{t) in ^ such that the field 
keeps moving to the right. Although the position of the interface is transported in 
the same direction as the magnetic field, it is lagging behind w.r.t. the postion of 
the ground state interfaces (p{t) (which moves from site 10 to site 12). At the same 
time the width of the interface is growing bigger. Also the energy of ip{t) is lagging 
behind w.r.t. the periodicity of the spectrum of H{t) (see Figure |4]). The difference 
with the ground state energy is steadily increasing, and it is an interesting question 
whether this difference will eventually saturate. In this computation, the minimum 
norm of ipit) remains above 0.99 for the first period {t < v~^) and above 0.97 for the 
second period {t < 2v~^)), so despite the non-adiabatic transport, ^p{t) is still well 
represented in the DMRG basis constructed from the low-energy spectrum of H{t). 
The overlap with the ground state ip{t) drops to 0.80 during the first period, and 
further to 0.71 during the second period. 

We have tested other magnetic field perturbations such as a smooth field 

^(^) = f E(l + - yt)))\L(x - vt) si 

X 

where 




forxG [L/2-l,L/2 + l] 
otherwise 



For L a multiple of 4, this is again a single-site perturbation for t = nv~ \ n = 
0, 1,2,.. ., and interpolating smoothly in between. In this case, the non-differentiable 
curves in Figure |4] become smooth, but no other qualitative differences are observed. 

Let us summarize the two regimes in terms of the two parameters v and B. We 
are in the adiabatic regime if Cv <C 1. Since for small B, the adiabatic constant 
C is inverse proportional to B, we require w <C -B. For large B we know that C is 
proportional to B and thus we want that v <^ 1/B. The regime where the initial state 
is stationary is simply given by the condition that B <Siv. 
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Figure 4. S^-magnetization profiles of for 0<t<r = 2i' ^ (left) and 

energy of the ground state and of (right, top) and energy difference between 
i/)(t) and the ground state (right, bottom) (A = 2, B = 0.5, v = 0.1, 200 time 
steps, S'^-profiles plotted for every 20 time steps). 



5. Summary 

We have studied the propagation of a magnetic domain waU in the ferromagnetic 
Heisenberg chain due to a time dependent magnetic field. This situation is of 
immediate interest in submicrometcr magnetic wires [21 13] and in estabhshing 
ferromagnetic gates [T]. In the simplest case, the magnetic field is localized near 
the center of the domain wall (say at site xq) and then moved along the chain to site 
Xi. If its velocity is small then we are in the adiabatic regime and the domain wall is 
shifted from xq to xi thereby preserving its shape. We find that there is an optimal 
region of the strength of B-fields for which the domain wall mobility is highest and is 
reduced for large and small B-fields. 

On top of an analytical study using the standard Adiabatic Theorem and based 
upon properties of the Heisenberg model that have been proved in recent years, we 
have performed a numerical study. The latter is a variation of the recently developed 
time-dependent Density Matrix Renormalization Group method. This allows us to 
follow the Schrodinger time evolution of the domain wall with high accuracy in the 
adiabatic as well as the non-adiabatic regime, and discuss in detail the dependence on 
the parameters of our model such as the strength and speed of the magnetic field and 
the coupling constant in the Heisenberg interaction. 

Appendix A. DMRG algorithms 

The DMRG algorithms for ground state [5] as well as time-dependent [71 [HI [SI [ID] 
computations arc well known, see also |30| and [6j. However, degeneracy and non- 
translation invariance of the ground states of the (translation invariant) XXZ kink 
Hamiltonian require a few non-standard adaptations. To lift degeneracy, we restrict 
to a sector of fixed total 5'^-magnetization. For a chain of even length, the ground 
state with total magnetization zero has an interface centered in between the two middle 
sites. This state can be grown in the usual way by successively inserting two sites in 
the middle and targeting the zero magnetization state at each step. For a kink state 
with non-zero magnetization, we grow the system initially in the zero magnetization 
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sector. If eventually a total magnetization with interface close to the middle is needed, 
we perform the finite system convergence sweeps in the new sector, with the zero 
magnetization state as initial trial state. If an interface close to one of the edges is 
needed, we grow the system at the right moment by targeting a magnetization sector 
which increases or decreases with one each step (effectively inserting two up or two 
down spins in the middle, thus shifting the kink to the right or left). The DMRG 
enlarging process implicitly assumes a translation invariant Hamiltonian. However, 
if the breaking of translation invariance is due to a local perturbation as in our case 
(with external magnetic fields located on one or two sites), we can grow the system 
targeting the zero magnetization kink, and add the perturbation for the finite system 
convergence sweeps where translation invariance is no longer needed. The DMRG 
state thus obtained converges indeed very rapidly to the true ground state ip{0) of the 
perturbed Hamiltonian ([5|) at time t = 0. 

In the standard time-dependent DMRG, adaptation of the basis is done by 
computing and truncating new reduced density matrices for the evolved state ip{t). 
Here the situation is different because we need to compute both the time evolved state 
ip{t) and the true ground state ip{t) of H(t) (in the same DMRG basis, naturally). 
Therefore we use the ground state (or several low- lying states) of H{t) to adapt the 
Hilbert spaces. The details of our algorithm are as follows. The time period [0, t] is 
first divided in n-r time steps of length S ~ r/n-r by discretizing the time evolution, 
i.e., the evolved states ip{t) are defined by 

^(0) =^{Q), 

ijjinS) ^ c-'*'^("'') V((" - l)*^) , n 1, . . . , . 

Each of these time steps is further divided into ut smaller steps of length St = S/ut 
for the Trotter decomposition, so 

and e-'-^T^^""^) is expanded by 

such that each factor acts on two sites which are successively represented exactly 
in a DMRG sweeping process. Note that only the interaction /iL/2.L/2+i is time- 
dependent. We shall call the so approximated state tpxinS). Let us start at time 
where V'(O) = (^(0). To apply e""'^^*^ to it we need to have ip{0) written in a DMRG 
basis that represents the low energy states of H{S). To this end we use ip{Q) as an 
initial trial state and apply standard finite system DMRG sweeps targeting the ground 
state ip{S) of H{d), and update ip{0) along the way. At the end of the sweeps we have 
a new DMRG basis for H{S), its ground state ip{S), and "0(0) expressed in the new 
basis. This is the adaptive part of our algorithm. Next we can apply c""'^^^^ to the 
new representation of ipiO) using the Trotter sweeping process and obtain the evolved 
state iPt{S)- By construction iPt{5) and (p{6) are expressed in the same DMRG basis, 
and their overlap and other quantities can be computed. The subsequent time steps 
proceed in exactly the same manner. 

Although the motivation for adapting the DMRG bases in this way is inspired 
by the question to analyze the adiabatic approximation, we have found that it is also 
very convenient to compute the time-evolved state in the non-adiabatic regime. In 
this case, we need to target more low-energy states, but we do not need to increase the 
block dimension, unlike standard time-dependent DMRG which needs approximately 
twice the block dimension of ground state DMRG to achieve the same accuracy [9] . 
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Appendix B. Parameter values and software 

Our problem is suitable for DMRG with small block dimension as its entanglement 
entropy [31] goes to zero for large intervals around the support of the magnetic 
field perturbation (it is exactly zero for the non-perturbed Hamiltonian due to its 
frustration free property). The algorithm converges up to machine precision to the 
true, known energies of the ground state and first excited state with the number of 
kept states (block dimension) as low as 16 for a chain of length L ~ 20. We carried out 
3 system sweeps to obtain convergence of the ground states and of the DMRG bases 
at each time step. For given v, the final time t = 1/v was divided in time steps of 
length 1 which were further subdivided into ut = 100 Trotter steps. For the ground 
state and adiabatic transport computations, we targeted the 3 lowest energy states of 
H{t) to construct the DMRG basis, increasing to the 5 lowest energy states for the fast 
changing magnetic field. A straightforward error analysis shows that |Ct — C| = C'(|)i 
where Ct it the adiabatic constant computed using the Trotter approximation to the 
time-evolved state, C is the true adiabatic constant, and ^ = l/{vn^)max{l,B^). 

We have implemented the ground state and time-dependent DMRG algorithms in 
a Matlab Toolbox which can be easily applied to other models as well. The software is 
included in the tar archive with the I^T[t;X source and figure files of the paper, available 
for download at arXiv : cond-mat/0702059. 
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